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Rotation and shear flows are ubiquitous features of many astrophysical and geophysical bodies. 
To understand their origin and effect on turbulent transport in these systems, we consider a forced 
turbulence and investigate the combined effect of rotation and shear flow on the turbulence proper- 
ties. Specifically, we study how rotation and flow shear influence the generation of shear flow (e.g. 
the direction of energy cascade), turbulence level, transport of particles and momentum, and the 
anisotropy in these quantities. In all the cases considered, turbulence amplitude is always quenched 
due to strong shear (£ = vky/A <C 1, where A is the shearing rate, v is the molecular viscosity and 
k y is a characteristic wave-number of small-scale turbulence), with stronger reduction in the direc- 
tion of the shear than those in the perpendicular directions. Specifically, in the large rotation limit 
(fl 2> A), they scale as A~ and -4 _1 | hi£|, respectively, while in the weak rotation limit (CI <gi A), 
they scale as A -1 and A~ 2 ^ 3 , respectively. Thus, flow shear always leads to weak turbulence with 
an effectively stronger turbulence in the plane perpendicular to shear than in the shear direction, 
regardless of rotation rate. The anisotropy in turbulence amplitude is however weaker by a factor 
of £ 1,/3 | ln£| (oc ^4 -1/,3 | ln£|) in the rapid rotation limit (Q. 3> A) than that in weak rotation limit 
(fl <C A) since rotation favors almost-isotropic turbulence. Compared to turbulence amplitude, 
particle transport is found to crucially depend on whether rotation is stronger or weaker than flow 
shear. When rotation is stronger than flow shear (fi A), the transport is inhibited by inertial 
waves, being quenched inversely proportional to the rotation rate (i.e. oc SI -1 ) while in the opposite 
case, it is reduced by shearing as A" 1 ■ Furthermore, the anisotropy is found to be very weak in the 
strong rotation limit (by a factor of 2) while significant in the strong shear limit. The turbulent vis- 
cosity is found to be negative with inverse cascade of energy as long as rotation is sufficiently strong 
compared to flow shear (fl S> A) while positive in the opposite limit of weak rotation (fi <C *4). 
Even if the eddy viscosity is negative for strong rotation (Q ^> A), flow shear, which transfers energy 
to small scales, has an interesting effect by slowing down the rate of inverse cascade with the value 
of negative eddy viscosity decreasing as \vt\ oc A~ 2 for strong shear. Furthermore, the interaction 
between the shear and the rotation is shown to give rise to a novel non-diffusive flux of angular 
momentum (A-effect), even in the absence of external sources of anisotropy. This effect provides a 
mechanism for the existence of shearing structures in astrophysical and geophysical systems. 

PACS numbers: 47.27. Jv,47.27.T-,97.10.Kc 

I. INTRODUCTION 

Rotating turbulent flows can be found in many areas such as engineering (turbo- machinery, combustion engine), 
geophysics (oceans, Earth's atmosphere) or astrophysics (gaseous planets, galactic and accretion disks). Large-scale 
fluid motions tend to appear as a robust feature in these systems, often in the form of shear flows (such as circulations 
on the surface of planets, differential rotation in stars and galaxies or flows in a rotating machinery). There have been 
accumulating evidence that large-scale shear flows as well as rotation play a crucial role in determining turbulence 
properties and transport, such as energy transfer or mixing (see below for more details). The understanding of the 
physical mechanism for the generation of large-scale shear flows and the complex interaction among rotation, shear 
flows and turbulence thus lies at the heart of the predictive theory of turbulent transport in many systems. 

A. Summary of previous works 

While both rotation and shear flow apparently have a similar effect on quenching turbulent transport, the efficiency 
of their effects as well as the basic physical mechanisms are totally different. It is thus useful to contrast these in 
detail. 

1. Sheared turbulence 

The main effect of shear flow is to advect turbulent eddies differentially, elongating and distorting their shapes, 
thereby rapidly generating small scales which are ultimately disrupted by molecular dissipation on small scales (see 
Fig. [1]). That is, flow shear facilitates the cascade of various quantities such as energy or mean square scalar density 
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FIG. 1: Sketch of the effect of shear on a turbulent eddy. 



to small scales (i.e. direct cascade) in the system, enhancing their dissipation rate. As a result, turbulence level 
as well as turbulent transport of these quantities can be significantly reduced compared to the case without shear. 
Another important consequence of shearing is to induce anisotropic transport and turbulent level since flow shear 
directly influences the component parallel to itself (i.e. x component in Fig. 1) via elongation while only indirectly 
the other two components (i.e. y and z components in Fig. 1) through enhanced dissipation. This shearing effect 
of shear flow can be captured by time-dependent Fourier transform where the wave number in the shearing direction 
(e.g. k x in Fig. I) increases linearly in time [l], 0, Hi- 
lt is important to emphasize that the aforementioned shearing effect (due to differential advection) is via nonlocal 
interaction between large and small-scale modes, and can dominate over nonlinear local interaction between small 
scales for sufficiently strong flow shear [e.g. H|. Therefore, the evolution of small-scale quantities can be treated as 
linear by neglecting local interactions compared to nonlocal interactions. This formulation, also called the rapid 
distortion theory (RDT) by various previous authors 0, HI, was used to study the linear response of turbulence to a 
mean flow with spatially uniform gradients. The linear treatment of fluctuations by incorporating strong flow shear 
was also used in the astrophysical context by [l[ by using shearing coordinates. The generation of large-scale shear 
flows (the so-called zonal flows) through a similar nonlocal interaction has been intensely studied in the magnetically 
confined plasmas, where turbulence quenching by shear flow is believed to be one of the most promising mechanisms 
for improving plasma confinement (fl 0] ■ 

In decaying sheared turbulence, |8( have shown a surprisingly good agreement between the RDT predictions and 
numerical simulations. Forced sheared turbulence was proposed for the first time by Q in the context of two- 
dimensional near- wall turbulence to explain the logarithmic dependence of the large scale velocity on the distance to 
the wall. In that case, the external forcing is provided by a continuous supply of vorticity from intermittent coherent 
burst of vorticity coming from the viscous layer. This work was later generalized to three dimensions (Tol . [ll| with 
the same conclusions. Subsequently, theoretical predictions (using a quasi-linear theory) for the transport of passive 
scalar fields in 2D hydrodynamic turbulence by [l2| and fl3| have been beautifully confirmed by recent numerical 
simulations In particular, they have shown that turbulent transport of particles can be severely quenched inversely 
proportional to flow shear A while turbulence level is reduced as A~ 5 ^ 3 . Ref. @ has shown that in 3D forced HD 
turbulence, strong flow shear can quench turbulence level and transport of particles with strong anisotropy (much 
weaker along the flow shear which is directly affected by shearing) and has emphasized the difference in turbulence 
level and transport, which is often used interchangeably in literature. A similar weak anisotropic transport was shown 
for momentum transport by [15| in forced 3D HD turbulence. Further investigations have been performed on turbulent 
transport in forced turbulence by incorporatin g th e interaction of sheared turbulence with different types of waves 
that can be excited due to magnetic fields [la Tl7l fl8| , stratification [li| or both magnetic fields and stratification 
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2. Rotating turbulence 

Rotation has both similar and different effects on turbulent transport. First, rotation can reduce transport in the 
limit of rapid rotation (similarly to flow shear), but through a physical mechanism that is different from that of shear, 
namely by phase mixing of inertial waves [2l| . It also induces only slight anisotropy in the transport (by a factor 
of two), much less significant than the strong anisotropy due to shear. Further, since phase mixing affects turbulent 
transport without necessarily quenching turbulence level, turbulence level may not be affected by rotation. This 
reduction in transport without much effect on turbulence level is a common feature of turbulence strongly affected by 
waves, and is also found in MHD turbulence where magnetic fields support Alfven waves Qj], [HJ and stratified 
turbulence [l9j where stable stratification excites internal gravity waves. A more striking difference between flow shear 
and rotation is that rotation facilitates the cascade energy to large scale, generating large-scale flows. For instance, in 
the extreme limit of very rapid rotation, the fluid motion becomes independent of the coordinate along the rotation 
axis (the so-called Taylor-Proudman theorem [23|,[24j]). The generation of large-scale flow has been shown by various 
numerical simulations including [25] and [26| . In particular, [26[ have shown that the inverse cascade of energy is 
more pronounced in forced turbulence due to statistical triadic transfer through nonlocal interaction. 

It is important to note that this nonlocal interaction leading to inverse cascade can be successfully captured by 
inhomogeneous RDT theory which permits the feedback of the nonlinear local interaction between small scales onto 
the large scales via Reynolds stress (constituting the other part of quasi-linear analysis) while neglecting nonlinear 
local interaction between small scales for fluctuations compared to nonlocal interactions. As must be obvious by 
comparing the Coriolis force with nonlinear advection terms, the RDT works well for sufficiently strong rotation 
(small Rossby number) even in the absence of shear flow. For instance, the agreement of the RDT prediction with 
numerical results has been shown by various previous authors including [25| . but mostly in decaying turbulence. 
However, in this case, the RDT cannot accurately capture the turbulence structure in the plane perpendicular to 
rotation axis where nonlinear local interactions between inertial waves seem important (see, e.g. [26]). The validity 
and weakness of the RDT together with comparison with various numerical simulation (without an external forcing) 
with/without shear flows and stratification can be found in excellent review by [13] and Cambon and [28], to which 
readers are referred for more details. 

In comparison, far much less is understood in the case of forced turbulence. In particular, the main interest in 
forced turbulence is a long-term time behavior where the dissipation, enhanced by shear distortion, is balanced by 
energy input, thereby playing a crucial role in leading to a steady equilibrium state. The computational study of 
this long time behavior is however not only expensive but also difficult because of the limit on numerical accuracy, 
as noted by (29| . Therefore, analytical theory by capturing shearing effect (such as quasi-linear theory with time- 
dependent wavenumber) would be extremely useful in obtaining physical insights into the problem as well as guiding 
future computational investigations. We note that the previous works by Kichatinov and Rudiger and collaborators 
using quasi-linear theory are valid only in the limit of weak shear. We further note that physically, 
the local nonlinear interactions in Navier-Stokes equation can be captured by an external forcing [35j, l36j . 

B. Main objectives and methodology 

Our main motivation is to understand the origin of large-scale shear flow and its effect on turbulent transport in 
rotating systems. To this end, we consider a forced turbulence and investigate the combined effect of rotation and 
shear flow on the turbulence properties including transport of momentum and particles. Specifically, we are interested 
in how rotation and flow shear influence the generation of shear flow (e.g. the direction of energy cascade), turbulence 
level, transport of particles and momentum, and the anisotropy in these quantities. Given the differences/similarities 
in the effects of flow shear and rotation (as discussed in Sec. II Ap . of particular interest is to identify the relative 
strength of flow shear to rotation rate for the cross-over between inverse and direct cascades and isotropic and almost- 
isotropic turbulence/transport. Recalling that flow shear of strength A acts over the time-scale A^ 1 while rotation 
induces inertial waves of frequency ~ Q, one could naively think that flow shear would dominate the effect of rotation 
for sufficiently strong shear with A 3> ^ while the effect of flow shear may be neglected in the opposite limit A <C f2. 
This will however be shown to be true only in the case of the transport of passive scalar fields and for the sign of 
eddy viscosity. That is, even in the case of weak shear compared to rotation A -^i fl, the shear has yet a crucial 
effect on determining the overall amplitude of turbulence level and momentum transport since its shearing process 
(generating small scales) works coherently over more than one oscillation of the waves. To complement this, we are 
also interested in how shear-dominated turbulence is influenced by rotation. As will be shown later, when the system 
is linearly stable, weak rotation tends to make turbulence/transport more 'isotropic'. 

Concerning momentum transport, another important question is the possibility of non-diffusive transport. In 
rotating turbulence, the inverse cascade can occur not only due to a (diffusive) negative viscosity, but also due to 
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non-diffusive momentum transport. The latter is known as the anisotropic kinetic a-effect (AKA) [37| or as the 
A-effect in the astrophysical community. The appearance of non-diffusive term in the transport of angular momentum 
prevents a solid body rotation from being a solution of the Reynolds equation [H, [3§| , and thus act as a source for the 
generation of large-scale shear flows. For instance, this effect has been advocated as a robust mechanism to explain 
the differential rotation in the solar convective zone. Starting from Navier-Stokes equation, it is possible to show 
that these fluxes arise when there is a cause of anisotropy in the system, either due to an anisotropic background 
turbulence (see (33j and references therein) or else due to inhomogeneities such as an underlying stratification. We 
will show that non trivial A-effect can result from an anisotropy induced by shear flow on the turbulence even when 
the driving force is isotropic, in contrast to the case without shear flow where this effect exists only for anisotropic 
forcing [32j |. 

We note that although much less attention has been paid to the effect of rotation and shear on mixing and transport 
of scalars (such as pollutants, heat or reacting species) compared to momentum transport, this is an important problem 
in understanding the distribution and mixing of a variety of physical quantities in different systems. For instance, 
observations show that the concentration of light elements at the surface of the Sun is smaller than what is expected 
by comparison with Earth's or meteorites abundance. As these light elements can only be destroyed below a strong 
shear layer (the so-called solar tachocline), their transport is subject to the effects of strong shear and rotation. The 
study of transport of passive scalar has been mostly limited to the purely rotating case [I(J, El[ or non- rotating sheared 
turbulence [43 . |43| ] . For purely rotating turbulence, linear theory has shown a strong suppression of particle diffusion 
by rotation, confirmed by numerical simulations [4l| . In comparison, the study of particle diffusion in sheared rotating 
turbulence was done only by [44| , who found that numerical simulation results agree fairly well with his linear theory. 

The purpose of this paper is to provide theoretical prediction on these issues by considering a 3D incompressible 
fluid, forced by a small-scale external forcing. As we are interested in the effect of flow shear, we capture this effect 
non-perturbatively by using time-dependent wavenumber [see Eq. (|3"])]. By assuming either sufficiently strong shear 
or rotation rate, we employ a quasi-linear analysis to compute turbulence level, eddy viscosity, and particle transport 
for temporally short-correlated, homogeneous forcing. As the computation of these quantities involve too complex 
integrals to be analytically tractable, they are analytically computed by assuming an ordering in time scales. In our 
problem, there are three important (inverse) time-scales: the shearing rate A, the rotation rate Q and the diffusion 
rate T> = i/ky where v is the (molecular) viscosity of the fluid and k~ x is a characteristic small scale of the system. 
We first distinguish the two cases of strong rotation (fl ^> A) and weak rotation (fi <C -4). The first regime of strong 
rotation will be studied in the strong shear (A S> T>) and weak shear (A <§; V) regime. On the other hand, the second 
regime of weak rotation will be considered only in the strong shear (A 3> T>) case, as the effects of both shear and 
rotation disappear in the opposite limit (A <^T>). We believe that our results would provide not only useful physical 
insights in understanding the complex dynamics of rotating sheared turbulence, but also serve as a guide for further 
theoretical/computational works, especially considering the difficulty of numerical study of this system. 

The remainder of the paper is organized as follows: in ^TTl we formulate our problem. Theoretical results of 
turbulent intensity and turbulent transport are provided in mill Some of the detailed analysis are provided only 
in i jllll We then discuss our findings in the strong shear limit in QYV\ and provide concluding remarks in <|V] The 
effect of rotation on linear stability of shear flows and some of the detailed algebra are provided in Appendices. Since 
analytical analysis performed in the paper are quite involved, some of the readers who are mainly interested in the 
results might wish to go to £HVI and <|V] after reading SjTTJ 



II. MODEL 

We consider an incompressible fluid in a rotating frame with average rotation rate f2, which are governed by 

<9iU + u-Vu = -VP + iA7 2 u + F - 2£l x u , (1) 
Vu = 0. 

Following 0], we study the effect of a large-scale shear Uo = Uo(x)j on the transport properties of turbulence by 
writing the velocity as a sum of a shear (chosen in the x-direction) and fluctuations: u = Uo + v = Uo(x)j + v = 
—xAj + v. Without loss of generality, we assume A > 0. In the following, we consider the configuration of Figure © 
where the shear and rotation (in the z direction) are perpendicular and simplify notation by using f2 = 2f2. Then, the 
Coriolis force is simply Q[— u y i + u x j], where i, j and k are the unit vectors associated with the Cartesian coordinates. 
Note that our x — y coordinates are not conventional in that our x and y directions correspond to y and x in previous 
works (see [29| for instance). Therefore, the shearing, the stream- wise and the span- wise direction correspond to the 
x, y and z direction, respectively. 
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FIG. 2: Sketch of the configuration in the perpendicular case 



To calculate turbulence amplitude (or kinetic energy) and turbulent transport, we need to solve the equation for the 
fluctuating velocity field. To this end, we employ the quasi-linear theory [45[ where the nonlinear local interactions 
between small scales are neglected compared to nonlocal interactions between large and small scales and obtain: 



d t v + U • Vv + v • VU 
V -v 



-Vp + ^V 2 v + f - ft x v 



(2) 







where p and f are respectively the small-scale components of the pressure and forcing. As noted in the introduction, 
this approximation, also known as the RDT is justified in the case of strong shear as the latter induces a weak 
turbulence, leading to weak interaction between small scales which is negligible compared to the (non-local) interaction 
between the shear and small scales. This has in fact been confirmed by direct numerical simulations, proving the 
validity of the predictions of quasi-linear theory with a constant-rate shear both in the non-rotating Q and rotating 
unforced [2!| turbulence and also for forced turbulence [l4j]- Further, note that the quasi-linear analysis is also valid 
in the limit of rapid rotation [3(| . 

To solve Eq. @ , we introduce a Fourier transform with a wave number in the x direction evolving in time in order 
to incorporate non-perturbatively the effect of the advection by the mean shear flow [H, 0, 9 : 



v(x,i) 



(2tt) 



d 3 k e i[k x {t)x+k y y+k zZ [-^^ ^ 



(3) 



where k x (t) 
velocity: 



k x (0) + k y At. From Eqs. Q and Q, we obtain the following set of equations for the fluctuating 



Ad T v x 
Ad T v y — Av x 
Ad T v z 




— ikyTp + f X + flVy 

-ikyp + f y - Vtv x , 
-ik z p + f z , 
tv x + v y + (3i) z . 



(4) 



Here, the new variables v = v exp[i/(kjjt + k x /3k y A)) and similarly for f and p have been used to absorb the diffusive 
term, and the time variable has been changed to r = k x (t)/k y . In the remainder of the paper, we solve Eq. ((4|) for 
the fluctuating velocity (with a vanishing velocity as initial condition) . We then use these results and the correlation 
of the forcing (defined in mi C|) to compute the turbulence intensity and transport (defined in mi Bp . 



A. Transport of angular momentum 



As the large-scale velocity is in the y direction, we are mostly interested in the transport in that direction. The 
large-scale equation for the y component of velocity Up is given by Eq. |T]) with a supplementary term V • R where 
R is the Reynolds stress given by: 

R = (v % ) . (5) 

To understand the effect of R on the transport of angular momentum, one can formally Taylor expand it with 
respect to the gradient of the large-scale flow: 

Ri = AiU - vrdjltfa + ■■■ = K- L U + v T A5n + ... . (6) 

Here, Aj and vt are the two turbulent transport coefficients from non-diffusive and diffusive momentum flux, re- 
spectively. Note that the first term in the expansion is due to the small-scale driving and the Coriolis force in Eq. 
(fTJ) which break the Galilean invariance 46]. First, is the turbulent (eddy) viscosity, which simply changes the 
viscosity from the molecular value v to the effective value v + vt- Note that the sign of eddy viscosity represents the 
direction of energy cascade, with positive (negative) value for direct (inverse) cascade. Second, the first term involving 
Aj in equation ([6]) is proportional to the rotation rate rather than the velocity gradient. This means that it does not 
vanish for a constant velocity field and thus permits the creation of gradient in the large-scale velocity field. This 
term bears some similarity with the a effect in dynamo theory [4?], [H| and has been known as the A-effect [3(J HI] 
or anisotropic kinetic alpha (AKA)-effect [13] • Similarly to the a effect, this effect exists only if the small-scale flow 
lacks parity invariance (going from right-handed to left handed coordinates). However, in contrast to the a effect, the 
A effect requires anisotropy for its existence [30, [37| • 



B. Particle (or heat) transport 

To study the influence of rotation and shear on the particle and heat transport, we have to supplement equation (fT|) 
with an advection-diffusion equation for these quantities. We here focus on the transport of particles since a similar 
result also holds for the heat transport. The density of particles iV(x, t) is governed by the following equation: 

d t N + U ■ VA = DV 2 N , (7) 

where D is the molecular diffusivity of particle. Note that, in the case of heat equation, D should be replaced by 
the molecular heat conductivity X- Writing the density as the sum of a large-scale component Nq and small-scale 
fluctuations n (N = Nq + n), we can express the evolution of the transport of chemicals on large scales by: 

d t N + U ■ VA = {p6n + B% )d i d j N , (8) 

where the turbulent diffusivity is defined as (i»jn) = —DjldjNa. D 1 ^ will analytically be computed to see the effect of 
rotation and flow shear on turbulent transport of chemicals which can be highly anisotropic. Note that the transport 
of a passive scalar quantity (contrary to the angular momentum which is a vector quantity) has to be diffusive due 
to the fact that it is solely advected by the flow [49j . 

For simplicity, we assume a unit Prandtl number D = v and apply the transformation introduced in equation ([3]) 
to the density fluctuation n to obtain the following equation: 

Equation ([9]) simply shows that the fluctuating density of particles can be obtained by integrating the fluctuating 
velocity in time. 



C. External forcing 



As mentioned in introduction, we consider a turbulence driven by an external forcing f . To calculate the turbulence 
amplitude and transport defined in Qll Al and ^11 Bl f which involve quadratic functions of velocity and/or density), we 
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prescribe this forcing to be short correlated in time (modeled by a ^-function) and homogeneous in space with power 
spectrum ip^j in the Fourier space. Specifically, we assume: 

(/iOci.tij/^ka.ta)) =T f (2ir) 3 6{k 1 + k 2 ) <f(t x - 1 2 ) ^-(k 2 ) , (10) 

for i and j = 1, 2 or 3. The angular brackets stand for an average over realizations of the forcing, and r/ is the (short) 
correlation time of the forcing. Note that the 6 correlation is valid as long as the correction time tj is the shortest 
time-scale in the system [i.e. r/ <§C fi" 1 ,^ -1 , l/(uk 2 )]. 

For most results that will be derived later, we assume an incompressible and isotropic forcing where the spectrum 
of the forcing is given by: 

^W^Fik^-hkj/k 2 ) . (11) 
It is easy to check that in the absence of rotation and shear, this forcing leads to an isotropic turbulence with intensity: 

2r f f°°F(k) 



dk , (12) 



where the subscript stands for a turbulence without shear and rotation. 

In addition to an isotropic forcing, we will also consider an anisotropic forcing in ^III A 21 to examine the combined 
effect of rotation and anisotropy, which can lead to non-diffusive fluxes of angular momentum. Specifically, we consider 
an extremely anisotropic forcing with motion restricted to a plane perpendicular to a given direction g. The motion 
in this perpendicular plane is however assumed to be isotropic. Such a forcing can be modeled by the following power 
spectrum [33j: 



rP ij (k) = G(k) 



kikj (g ■ k) 2 g-k 
Oij - -p ^— dij - g l g j + -p-(5iKj + 9jh) 



(13) 



In that case, the turbulence without rotation or shear would have the following properties: 

((vo-g) 2 ) - 0, (14) 
// ^2\ 2r / [°° G ( k ) ,1 

((voxg)) = 3wX — dk - 



III. ANALYTICAL RESULTS 



The system (j4} can be simplified to: 

d T v z 



M - 



hi(r) 



A 



■-dr[TV x \ +13 V x 

7 7 
-{tv x + (3v z ) . 



-0/3 



A 



(15) 



Here: 



n = 

hi(r) = 



fl/A , (3 = k z /k y 
1 fx - rf y - Prf z , 



7 = 1 + P 2 = k 2 H /k 
h 2 (r)=f z -Pf y . 



(kjj — k 2 + k z ) , 



(16) 



To solve the first of equation (j 1 5[) which is a non-homogeneous second order differential equation, we need two 
boundary conditions. We impose a vanishing initial velocity v(to) = which implies v x (tq) = and d T v x \ T=ro = 
^i(t"o)/(7 + Tq )A. The second boundary condition can be shown to be obtained in the intermediate steps of deriving 
Eq. G5j. 

The exact solution to (fTSf is obtained in the appendix|Al where we address the stability of the homogeneous solution 
of system (with f = 0). Computations of correlation functions, by using this exact solution, however turns out to be 
too complex to be analytically tractable. To gain a physical insight into the role of inertial waves and flow shear in 
turbulent transport, we consider the two limits - (i) the strong rotation where the effect of waves dominates shearing 
(O 3> ^4) and (ii) the weak rotation where shearing dominates the effects of waves (S7 <C -4). Approximate solutions 
can be derived in these two regimes which can then be used for deriving analytic form of correlation functions for 
turbulence intensity and transport. 
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A. Rapid rotation limit: Q ^> A 

When the rotation rate is much larger than the shearing rate (fio = |f2|/-4 3> 1), the oscillation of inertial waves is 
roughly coherent without being damped over shearing time of A~ 1 . Therefore, these waves can play a dominant role 
in determining the direction of energy cascade (sign of eddy viscosity) and transport of particles via phase mixing 
(i.e. by affecting the phase relation). However, as shown below, flow shear can still have a non-trivial effect on 
turbulence by enhanced dissipation so long as it is stronger than molecular dissipation. To characterize the latter, 
we introduce a parameter £ = vk 2 /A, the ratio of typical molecular dissipation rate to shearing rate. Here, k y is 
the characteristic wavenumber of the forcing in the stream-wise direction. We can, for instance, envision the forcing 
to have a spectrum peaked around this characteristic wave-number k y . In the following, we examine the changes in 
turbulence characteristics in weak (£ ^> 1) and strong ((< 1) shear limits to elucidate the effects of flow shear in 
inertial wave-dominated turbulence. 

In the rapid rotation limit (|f2| ^> A), the solution of Eq. (fTS)) can be found by using WKB approximation as: 

Mr) = A{i + \ 2)3/4 jf dt | (7 h + f )1/4 cobW, r)] + k(t)(7 + t 2 )^9sin[v(t, r)]| , 

^ (T) = A^ + r^)3/4 £ dt { ( 7 + ^i/4 (~ T cos ^ t)1 + WVTT^smMt, r)]) 

+h 2 (t)(j + t 2 ) 1/4 {-9t sin[u(t, r)] - /3^7 + r 2 cos[«(t, r)]) } , (17) 

%{r) = Al{l l T2f/i jf dt{ (<y H + f )1/4 (-fh co S [v(t, r)] - e^T^sm[v(t, r)]) 
+ha(t)(rf + i 2 ) 1/4 (-O0T ain[v(t, t)] + ^7 + t 2 cos[v(t, r)]) } . 



Here, 



ft = |fi| , lo = \0\SI O , 6 = siga(j30,) , (18) 

I arrsinVi ( — — I J- D f — 



s(t) = 1 arcsinh + O . 



v(t,r) = iv [a(t) - s(r)] . 

In the following subsections, we compute the various correlation functions by assuming a homogeneous and short- 
correlated forcing [see Eq. (TlTJ|) ]. As the system (fT5|) involves the forcing in terms of hi and h 2 only [see Eq. (fTS)) ], it 
is convenient to use the power spectrum <fiij as: 

(^(ki,ii)^(k 2 ,t 2 )) =r / (2^) 3 «5(k 1 +k 2 ) ( 5(t 1 -t 2 )^(k 2 ) , (19) 

for i and j = 1 or 2. In the case of an isotropic and incompressible forcing [Eq. (|11[) ]. in Eq. (|19p can be written: 

0n(k) = 7(7 + a 2 )F{k) , 12 (k) = , 22 (k) = jF(k) . (20) 



_Z. Turbulence intensity 



We begin by examining the effects of rotation and flow shear on turbulence level in wave-dominated turbulence 
due to strong rotation (|f2| ^> .4). The effect of shear will further be clarified by comparing results in weak shear 
limit (£ 3> 1) with those in the strong shear limit (£ <C 1). First, turbulence intensity in the shear direction can be 
obtained by using Eqs (fT7]) and (fT9j) as: 



>12 



(k) sin[2u(o,r)] + </> 22 (k) + a 2 sin 2 [v (a, t)] | 
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Here, a — k x /k y , (3 — k z /k y , 7=1 + /3 2 , £ = (yky)/A and Q(x) = x 3 /3 + jx. In the case of an isotropic forcing [Eq. 
(|2H)) ] , Eq. (f2"Tj) and the turbulence intensity in the two other directions can then be derived as: 



(vl) Tf 



(2tt)M 



d 3 k 1 ^ 1 + a 2 F{k) 2°(k) , (22) 



Here: 



(«2> = (^4/^^/7+^^) {/3 2 /°k)+/ 2 (k)} , 
(^ 2 > = (^4 / d3fc {/°(k)+/ 2 (k)} . 



r+oo p -2£[Q(r)-Q(a)] 



(23) 



In order to understand the effect of shearing on turbulence intensity in this wave-dominated turbulence, we first 
examine (|22[) in the weak shear limit (£ 3> 1) where the shear is negligible. In this case, the integral I p in Eq. (|23|) 
takes the approximate value: 

p (It) = — (24) 

V ' 2£( 7 + a 2 ) 5 / 2 2^fc 2 ( 7 + a 2 ) 3 / 2 ' V ; 

By using Eq. (124)) in Eq. {22}, we can then obtain the following result for the turbulent intensity: 

<«2> = 7 ^[*k?3-±- t (25) 



to* 



(2?r) 3 J 2vk 2 1 + a 2 



y/ (2tt) 3 J 2vk 2 7 + a 2 ' 



v 2 ) = / d * k 



,2 



(2tt) 3 / 2^fc 2 7 + a 2 ■ 
Performing the integration over the angular variable, we obtain: 

(vl) = 7 ll— [ dk^^- [ d(f> [ d9 sin 9 (cos 2 + sin 2 6 sin 2 



(27t) 3 7 2v J 







dk-^- = -(v 2 ) , (26) 



3(2tt) 2 J v 3 
1 



Here, (uq) is the turbulence amplitude in the absence of rotation and shear [sec Eq. ([f~2j) ]. These results thus show 
that, in the large rotation limit, the turbulence intensity is isotropic and equals to the one without rotation [see Eq. 
(|12p] for sufficiently weak shear with £ ^> 1. Furthermore, in this limit of a sufficiently weak shear where (SI, T>) ^> A, 
turbulence intensity is independent of rotation since waves do not necessarily quench turbulence level. A similar result 
was also obtained in MHD turbulence and stratified turbulence where magnetic fields and gravity waves mainly affect 
transport without much effect on turbulence level [13, H3| • We shall show below that a strong anisotropy can be 
induced when shearing effect is not negligible (£ -C 1) even in the rapid rotation limit (SI ^> .4). 

In order to understand the effect of flow shear, we now consider the strong shear limit (£ <C 1). In this limit, the 
integral (f2"3")) is simplified as: 

(27) 
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By plugging Eq. J27J) in Eq. ([221 we obtain: 



(vl) = J d*k V^F(k) ex t(v 2 o) , (28) 



y/ (2tt)M 



(vl) = (2^4 / ^ v/t+^W cx £|ln£|(z; 2 }, 

to leading order in £ <C 1. Note that in the calculation of (v 2 ), we neglected the component proportional to a = k x /k 
as it is odd in both k x and k y and thus vanishes after integration over the angular variables for an isotropic forcing. 
The last terms in Eq. (|28p. expressed in terms of the turbulence amplitude in the absence of rotation and shear (vq) 
[see Eq. (TT21) ]. explicitly show the dependence of turbulence level on rotation and shear. That is, all the components 
of turbulence intensity is reduced for strong shear £ <C 1. Further, the x component along shear is reduced as £ oc A^ 1 
while the other two components as £| ln£|, with an effectively weaker turbulence in the shear direction than in the 
perpendicular one, by a factor of ln£. This shows that shear flow can induce anisotropic turbulence (unlike rotation) 
even when the forcing is isotropic. This result is similar to that obtained in the simulation of a Couette flow at high 
rotation rate [5l| where the velocity fluctuations perpendicular to the wall exceed that in the stream-wise direction. 
Nevertheless, Eq. (f2"B")) shows that a strong rapid rotation yet insures an isotropy in velocity fluctuations in y — z 
directions {(v 2 ) — (v 2 z )). 



2. Transport of angular momentum 

As noted in the Introduction, rotation tends to cascade energy to large scales while shear flow to small scales. 
Would thus the inverse cascade be a robust feature for rapid rotation (fi ^> A)? If yes, what would the effect of 
flow shear? Would there be a non-diffusive momentum transport? We answer these questions by first considering an 
isotropic forcing and then anisotropic forcing. The effect of shear will be elucidated by looking at the two limits of 
weak shear (£ 3> 1) and strong shear (£ <C 1), as done in §3.1.1. 

First, in the case of an isotropic forcing, we obtain the following Reynolds stress from equations ([T7|) and (TT9l) : 



(v x v v ) - 



(2tt)M 



J d 3 k ^7 + a 2 F(k) / x (k) , (29) 



where I 1 was defined in Eq. ([23]) . Eq. (|29]) is computed in the weak and strong shear limits, below. 

First, in the weak shear limit (£ 3> 1), there is no contribution to leading order in fi -1 as the function I 1 is odd in 
a and thus vanishes after integration over the wave vector. We thus include one higher order in f2 _1 in the expansion 
and obtain the following result: 

<"-"»>-(^/ d3t ^ J(k >- <»> 

Here, we defined a function J(k), which has the following asymptotic behavior in the weak shear limit: 



J W = / (7 + t2) 3/2 sm[2u { S (a)-s(T)}]dT (31) 

aZJoA 



2(7 + a 2 ) 3 /2[i/ a *4+wg a 



where loq = ujqA/ y/ 7 + a 2 . Plugging Eq. (|31|) in Eq. (|30|) and performing the integration over the azimuthal angle 
variable ef>, we obtain: 

M = ^fAl I dk k2p ( k ) r de sin5 6 2m!u— 2 • ( 32 ) 
32(37r)^ J J v z k 4 +uj z 

Finally, we change the integration variable from 6 to ujq = £1 cos 9, obtaining the following formula: 



r A f +oc (l — W 2 IVt 2 Y 

^ = dkk2F ^l d ^ \^\ J • (33) 
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Therefore, in the large rotation and weak shear limit, the Reynolds stress becomes purely diffusive (with no A-effect) 
with the turbulent viscosity: 

32(27r) 2 |0| 7 v y ' 

This result shows that the turbulent viscosity is positive and proportional to fl" 1 for large fl. It is worth comparing 
Eq. ([34]) with Eq. (22) in [3l|. To this end, we use Eq. (fl"2|) . which gives the turbulence amplitude without rotation 
(the original turbulence of Kichatinov) in Eq. (|34]l to obtain the turbulent viscosity Vt ~ 7r(wo)/64|f2|. Thus vt in 
Eq. (|3"4")1 is the same as Eq. (22) in [3l| for |f2| 3> 1 and 6 = tt/2, but has an opposite sign. This is due to the 
r-approximation used by Kichatinov which gave an unphysical result. Later, [52j showed that the viscosity is also 
positive at any rotation rate when derived consistently with quasi-linear approximation in the weak shear limit. 
In comparison, in the strong shear limit (£ -C 1), the function I 1 in Eq. (f2"3"| has the following asymptotic behavior: 

I 1 M = —±=z. (35) 

Vi + a 

Plugging Eq. ()35|) in Eq. (|2"9")l . we obtain the turbulent viscosity in the strong shear limit as: 

Eq. (f36|) shows that the turbulent viscosity is negative (as F(k) > 0) in the strong shear limit, in sharp contrast 
to the weak shear limit where vt > [see Eq. (|34l) ]. Furthermore, the magnitude of i>t is reduced by the shear 
(oc A~ 2 ) and is independent of rotation, which should also be compared with the weak shear limit [see Eq. (f34|) where 
vt oc f2 -1 ]. Therefore, the turbulent viscosity changes from positive (for weak shear) to negative (for large shear) as 
the ratio of shear to dissipation increases. This result can be understood if we assume that, as in most rapidly rotating 
fluid, the inverse cascade is associated with the conservation of a potential vorticity [53j |. In the presence of strong 
shear (compared to dissipation), the potential vorticity is strictly conserved giving rise to an inverse cascade (negative 
viscosity). When the dissipation increases, the potential vorticity is less and less conserved and thus the inverse 
cascade is quenched. Our results show that there is a transition from inverse to direct cascade as the dissipation is 
increased. A similar behavior is also found in two-dimensional hydrodynamics (HD) where an inverse cascade can be 
shown to be present only for sufficient weak dissipation (l6| . 

It is important to note that the negative viscosity vt < obtained here for strong rotation/strong shear (fl 3> 
A vky) signifies the amplification of shear flow as the effect of rotation favoring inverse cascade dominates shearing 
(generating small scales). However, the magnitude of vt is reduced by shear as |i*r| oc A~ 2 since flow shear inhibits 
the inverse cascade. This can be viewed as 'self-regulation' - that is, self-amplification of shear flow is slowed down 
as the latter becomes stronger. 

The preceding results [Eqs. (f3"4")l and (|31)]) ] indicate that in the large rotation limit where rotation dominates over 
shear, the momentum transport is purely diffusive for isotropic forcing, with opposite sign of turbulent viscosity for 
weak (£ 3> 1) and strong shear (£ <C 1) for a fixed value of |f2|/.4 (3> 1). In the case of anisotropic forcing, there is 
however a possibility of the appearance of non-diffusive momentum transport (A-effect). To examine this possibility, 
we now consider an extremely anisotropic forcing (introduced in £|II CP where the forcing is restricted to horizontal 
plane (y-z), perpendicular to the direction of the shear. Using Eq. (flUl) with cjij = Sn, we obtain the following 
Reynolds stress: 

M = -T^O / d " k , lG[k \ [{ Jl ( k ) - J'(V}+f36K(k)] . (37) 
Here, I 1 was defined previously in Eq. (|23[) and: 

r+oc -2£[Q(T)-Q(a)] 

J '( k ) = J a ' (7 + r2)3 /2 cos [ 2w o - S ( T )>] dT . ( 38 ) 

r+oc _2£[Q(r)-Q(o)] 



K(k) = (7 + r2) sin [2c {s(a) - s(r)}] dr 



We again consider the weak and strong shear limits in the following. First, in the weak shear limit (£ ^> 1), Eq. 
37|) is simplified to: 

M - 7^7 / d " k J G l k) Z* • (39) 

y (2tt) 3 A J 4(7 + a 2 ) 3 / 2 v 2 k 4 + ujQ 
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Performing the angular integration in Eq. (J39J) and taking the large rotation limit, we obtain the following: 




3(27r) 3 04 



/ 




(40) 



Equation (|40|) is odd in the rotation and thus represents the A-effect. Again, the latter favors the creation of velocity 
gradient rather than smoothing it out and can thus provide a mechanism for the occurrence of differential rotation 
(e.g., in the sun). By using Eq. (fTl"]) . one can see that the A-effect is proportional to the anisotropy in the turbulence 
without shear and rotation. This result shows that, in the large rotation limit, one needs anisotropic forcing to generate 
non-diffusive fluxes of angular momentum [as in the case without shear as shown [3l[ . This should be contrasted to 
the case of weak rotation (see mil B[) where the shear can alone give rise to an anisotropic turbulence, thereby leading 
to a A-effect even with an isotropic forcing. 

Finally, in the opposite, strong shear limit (£ <C 1), Eq. ([57)) becomes: 



which is even in the rotation. Thus, the turbulent viscosity vt is obviously negative. Thus, in the large shear limit 
(but still negligible compared to the rotation), anisotropic forcing does not induce any non-diffusive fluxes but just 
increases the magnitude of the negative turbulent viscosity. 



In the large rotation limit (|f2|/^4 3> 1), inertial waves might play a crucial role in transport of particle as waves 
can alter the phase relation between particle density and velocity, as noted previously. How does this effect appear in 
forced turbulence? What is the effect of shear flow on particle transport dominated by waves? These questions are 
answered in this subsection. 

In the rapid rotation limit (|f2|/.4 3> 1), turbulent particle diffusivities can be obtained after a long, straightforward 
analysis (see Appendix [Bl for details about the algebra) as: 



Note that in that case, the result is not sensitive to the value of the parameter £ and thus we do not distinguish 
between the weak and large shear limits. Eq. I]42p shows that D^ x , and D" are all reduced as f2 _1 (with no 
effect of the shear) for large and also that there is only a slight anisotropy in the transport of scalar: the transport 
in the direction of the rotation is twice larger than the one in the perpendicular direction [34| . Interestingly, this 
anisotropy in the transport of particles is not present in turbulence intensity [see Eq. (|26|) ]. This is because waves 
can affect the phase between density fluctuation and velocity, not necessarily altering their amplitude. However, it 
is important to note that this anisotropy is only a factor of 2, much weaker than that in sheared turbulence without 
rotation Q. 

To summarize, in this subsection 3.1, we have examined how a shear flow can affect the turbulent property when 
turbulence is largely dominated by inertial waves in rapid rotation limit (|f2|/.4 1). In particular, the results show: 

1. that shear flow reduces turbulence level with a strong anisotropy [Eq. [)28p]. leading to an effectively weaker 
turbulence in the direction of the shear [which would otherwise be almost isotropic [Eq. (|25p]: 

2. that in comparison, transport of particles is mainly governed by waves with almost isotropic property (within a 
factor of 2) and quenched as fl" 1 as rotation rate fl increases; 

3. that energy cascade is inverse with negative viscosity for strong rotation/shear limit (tt ^s> A vk^) while its 
rate is slowed down by strong shear; 




(41) 



3. Transport of particles 




(42) 



4. that momentum transport is purely diffusive for isotropic forcing, with non diffusive transport appearing only 
for anisotropic forcing. 
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B. Weak rotation limit: f2 -C A 



When <C A, flow shear can distort inertial waves over the period of their oscillation, dramatically weakening the 
effects of these waves on turbulence. Therefore, shear may take a dominant role in determining turbulence property 
(studied in Q) while rotation modifies some of the properties of this shear- dominated turbulence. The investigation 
of this limit would thus permit us to clarify the effects of rotation as well as flow shear, thereby complementing the 
analysis done in Sec. IIII Al for strong rotation (£l^> A). Of particular interest is (1) to what extent the quenching and 
anisotropy of sheared turbulence 0] are affected by rotation, which favors isotropic turbulence; (2) how the direction 
of the energy cascade, which tends to be direct in 3D sheared turbulence, is affected by rotation (which prefers inverse 
cascade); (3) whether momentum transport can occur via non-diffusive fluxes. 

To answer these questions, we expand various physical quantities in powers of f2o = as: 

X(T)=X Q (T)+n X 1 (T) + ... , (43) 

in the weak rotation limit (f2 <C A) and calculate the turbulence intensity and transport up to first order in f2o- 
For the sake of brevity, we here just provide the final results of the calculation. Note that in this limit, we are only 
interested in strong shear case (£ <C 1) since in the opposite limit where vk^ ^ A^ Vt, the effects of both shear and 
rotation simply disappear to leading order. 



Turbulence intensity 



By using the expansion in powers of Hq (|43|) and Eq. (fT9| and after a long, but straightforward algebra, we can 
obtain the turbulence intensity in the shear direction as follows: 



(27r)M 



d 3 fc0 n (k) [L„(k)+/? 2 nLi(k)] 



(44) 



Here: 



r- + oo -2£[<3(r)-Q(a)] 

L (k) = / dT ; — dT 



( 7 + r 2)2 
+ oo -2f[Q(r)-Q(o)] 



( 7 + T 2 ) 



/•-l-oo 

£i(k) = / dr 

J a 

rr, N 1 ( x \ 

1 (X) = arctan 



r{T(T)-T(a)}-iln(l±4 
2 V 7 + a 1 



In the strong shear limit (£ <C 1), the integrals L and L x in Eq. (|43)l can be simplified: 



(45) 



dr , 



L (k) 

Li(k) 



+°° 1 _ _L_ 

(7 + r 2 ) 2 2^ 



7r , . a 

T(a) 7t 

2^7 7 + a 2 



+ OC 



1 



(7 + r 2 ) 2 



r{T(r)-T(a)}-hn(l±^ 
2 \7 + a / 



(46) 



f/r 



+ OC 



2 7 ( 7 ■ 



1 

2T 



T(r) 



{T(r) - T(a)}dr . 



Note that the second formula for Li in Eq. J46|) was obtained by integration by part. The leading order behavior of 
Eq. (|44p coming from the term involving L is due to shearing effect, showing that (u 2 ) is quenched by flow shear 
ex A^ 1 (see 3]). The effect of rotation appears as a correction proportional to L\. One can see from Eq. (|4"6")) that 
this correction L\ is positive for all values of a (for a < 0, the negative part of the integral is always smaller than the 
positive one as the first term is odd in r and the second one is an increasing function of a). Therefore, the turbulence 
intensity (v^.) in Eq. (|44| increases for > whereas it decreases for f2 < 0. This can physically be understood from 
the linear instability analysis (performed in appendix): that is, instability (fi > 0) increases turbulence level while 
stability (CI < 0) reduces it. 
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The other components of the turbulence amplitude can be obtained by following similar analysis in the strong shear 
limit (£ <C 1) as follows: 



(2tt)M 



(2tt)M 



d 6 k 



B 1 



T(a) 0n(k)+0 22 (k) 



2^ 

x [r(l/3) + n/3 2 r(4/3)(-lnO] 

2 



2^ 



T(o) 0n(k)+0 22 (k) 



3 7 2 U 



1 

3? V2^ 



1/3 



1/3 



(47) 



x [r(l/3) + n/3 2 r(4/3)(-ln£)] 



Here, T is the Gamma function. The first terms in Eq. (|47[) represent the turbulence amplitude in the direction 
perpendicular to shear without rotation [3], which are reduced as A~ 2 I Z for strong shear. Compared to the leading 
order behavior of (i> 2 ) oc A^ 1 in shear direction, the reduction is weaker by a factor of i; 1 / 3 . That is, a strong 
anisotropy in turbulence level can be induced for strong shear. The second terms in Eq. (|4T|) capture the effect of 
weak rotation on sheared turbulence, with turbulence amplitude again being increased or decreased depending on the 
sign of Cl. Furthermore, the correction comes with a multiplying factor oc |m£| > 1, which is larger compared to 
that for the amplitude in the shear (x) direction (which is independent of shear [Eq. (I44p ]~). Therefore, in the stable 
situation (fl < 0) of our interest, weak rotation has the effect of reducing turbulence in the y — z plane more than the 
one in the shear direction. As a result, the anisotropy induced by flow shear is weakened by rotation. Interestingly, 
this illustrates the tendency of rotation of leading to almost isotropic turbulence. 

It is also interesting to note that the leading order terms in (u 2 ) and (v 2 ), although apparently very similar, are 



not exactly the same. For instance, in the case of an isotropic forcing, the angular integration gives 



> (vi 



This 



slight anisotropy in y — z (stream and span-wise) directions in sheared turbulence was also observed in numerical 
simulations of homogeneous turbulence subject to high shear rate: the fluctuating velocity in the direction of the flow 
is larger than the one in the direction of the shear Q. This can be contrasted to the exact equipartition between (u 2 ) 
and (i> 2 ) [see (|28p] in the case of rapid rotation. This is another manifestation of the difference between shear flow 
and rotation in inducing anisotropic turbulence. 

In summary, in the case of a weak rotation/strong shear turbulence (A |fl| and A ^S> vk 2 ) : the rotation tends to 
reduce the anisotropy in sheared turbulence. 



2. Transport of angular momentum 



As noted previously, a strong anisotropy in turbulence is caused by strong shear in the weak rotation limit. There 
is thus a possibility that this anisotropic turbulence gives rise to non-trivial non-diffusive momentum transport. This 
will be shown to be the case below. 

In the strong shear limit (£ <C 1), momentum flux can be derived as: 



(k) 



1 



2( 7 + a 2 ) 



■0 



2 ( K 



\2^7 



T(a) 



(48) 



+ ^(-ln£) 
37 



( n 



V2V7 



T(a) <Mk)+0 22 (k) 



The momentum flux in (|48|) consists of a diffusive part (the first half term in the integrand on the RHS) and a 
non-diffusive part (the second half term in the integrand on the RHS). First, the diffusive part, independent of f2, 
recovers the eddy viscosity of sheared turbulence without rotation [3j, showing that its value decreases as oc A~ 2 for 
strong shear. This result agrees with previous studies of non-rotating sheared turbulence [Io| which found a Reynolds 
stress inversely proportional to the shear, leading to a log dependence on the distance to the wall for the large-scale 
shear flow. Second, the non-diffusive part, the correction due to the rotation, is proportional to Cl and is odd in the 
rotation. This is a non-diffusive contribution to Reynolds stress - the so-called A-effect. The origin of this non trivial 
A-effect is the strong anisotropy induced by shear flow on the turbulence even when the driving force is isotropic. It 
is important to contrast this to the case of rapid rotation limit where non-diffusive fluxes emerge only for anisotropic 
forcing. A similar result was also found in ijlll A 21 [see Eqs. (|3~4"|) and (|3T>)) ]. This A-effect [the second term in Eq. (|48[) ] 
is obviously of the same sign as 17 whereas the turbulent viscosity [the first term in Eq. (|48[) ] can either be positive 
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or negative, depending on the relative magnitude of the two terms inside the integral. In the two-dimensional (2D) 
limit with k z — ((3 — 0), we can easily show that the turbulent viscosity is negative. Note that in this 2D case, 
vt < signifies the amplification of shear flow while \vt\ oc A~ 2 reflects that the generation of shear flow slows down 
for strong shear. In contrast, in 3D with an isotropic forcing, the turbulent viscosity is positive. Finally, we note that 
our results here are compatible with previous studies which showed that non-diffusive fluxes of angular momentum 
[30L l3ll ] are proportional to the anisotropy in the background turbulence, which is induced by flow shear in our case. 



3. Transport of particles 

Transport of particles has been shown to be severely quenched by shear flow with strong anisotropic properties [3] • 
We now examine how (weak) rotation affects this. In the strong shear limit (£ <C 1), we can find turbulent diffusivity 
of particles as: 



D X T X 



D Z T Z 



(2n) 3 A 2 



d 3 k (£u(k) 



2^7 



T(a) 



i + n0 



(49) 



{2-KfA 2 



1 



2/3 



>22 



(k) 



HiJ r(2/3)|l + 2^- 

The first terms in Eq. (|49p manifest the quenching of particle transport for strong shear as DZj? oc A~ 2 and 
D?j? <x ^4 -4 / 3 , with effectively faster transport in span- wise direction compared to shear direction. That is, a strong 
anisotropic transport can arise for strong shear. It is interesting to contrast this result to that in the case of rotation 
where the transport in the shear (x) direction was larger only by a factor 2 than the one in the perpendicular direction. 
The second, correction terms in Eq. (|49]) represent the effect of rotation and are proportional to f2: Thus, for ft > 0, 
the transport is increased whereas it is reduced for O < 0. This is physically because a weak rotation destabilizes 
sheared turbulence for f2 > whereas it stabilizes for ft < (see figure 0] and the discussion in Appendix [A} . Note 
that a similar behavior was also found in turbulence intensity, given in Eqs. (|4"4"|) and (|4"7|) . Thus, one can see that 
for stable configuration (fi < 0) of our interest, the corrections due to rotation tend to weaken the strong anisotropy 
induced by flow shear. These results highlight the crucial role of shear in transport, in particular in introducing 
anisotropy. 

To summarize Sec. 3.2, in the slow rotation limit where turbulence is mainly governed by flow shear, turbulence 
intensity [Eqs. (|44)l and (|47|) ] and transport [Eq. ([49]) ] can be severely quenched with strong anisotropy due to 
shearing while weak rotation weakens this anisotropy to next order. The strong anisotropic turbulence was shown to 
give rise to a A-effect for momentum transport [Eq. (|48[) ] even for an isotropic forcing. 



IV. DISCUSSION 



In ^III[ depending on the values of the parameter £ = {yk 2 )/A, we considered two regimes: the strong shear (£ <C 1) 
and the weak shear limits (£ 3> 1). Since we are interested in the effects of flow shear as well as rotation, we here 
summarize and discuss our results obtained in the limit of strong shear with £ <C 1. Table [J summarizes our findings 
by highlighting the quenching of these quantities due to large shearing rate A and the rotation rate f2 (or their ratio, 

= Q/A). These results are discussed in the following. 



A. Turbulence amplitude 

In all the cases considered, turbulence amplitude is always quenched due to strong shear (£ = vk 2 /A < 1), with 
stronger reduction in the direction of the shear (x) than those in the perpendicular directions. Specifically, in the 
large rotation limit, they scale as A^ 1 and -4 _1 | ln£|, respectively while in the weak rotation limit, they scale as A^ 1 
and „4~ 2 / 3 , respectively. Thus, flow shear always leads to weak turbulence with an effectively stronger turbulence 
in the plane (y-z) than in the shear direction, regardless of rotation rate. The anisotropic reduction of turbulence 
amplitude is because of the shear which increases the dissipation (anisotropically) by efficiently creating small-scale 
fluctuations in the x-direction, with a direct impact on turbulence in the shear direction (see Figured]). The anisotropy 
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n > .4 


n < .4 




4- 1 


4- 1 [l + cn] 




•A -1 |ln£| 


A -2/ 3 [i + cn\ki£\] 




-.A" 2 


A' 2 


A* 





A~ 2 \lnt\ 






A~ 2 [l + Cn|ln£|] 


Df, y ~ D Z T Z 


n- 1 


a~ 4/3 [i + cn[ineO 



TABLE I: Summary of our results obtained in the strong shear limit (£ = vk y / A <C 1). The C symbol stands for an additional 
constant of order 1. 

in turbulence amplitude is however weaker by a factor of £ 1//3 | ln£| (oc ^l _1 / 3 | ln£|) in the rapid rotation limit than 
that in weak rotation limit since rotation favors almost-isotropic turbulence. In the case of weak rotation, the effect 
of shear on turbulence amplitude can be understood in terms of stability of rotating shear flow (see Appendix [5] for 
more details). In the case of weak rotation (f2 A), the effect of rotation appears in combination with the linear 
instability criterion in turbulence amplitude with linear stability f2 < (instability Cl > 0) decreasing (increasing) 
turbulence amplitude. For stable configuration S7 < 0, the rotation thus has the effect of weakening the anisotropy 
caused by strong shear. In summary turbulence amplitude is quenched by shear with strong anisotropy while rotation 
tends to weaken the shear-induced anisotropy. 

B. Transport of angular momentum 

The transport of angular momentum was found to involve two contributions: the turbulent viscosity vt and the 
A-effect. The former is a diffusive flux, making the effective viscosity to vt + v (y is the molecular viscosity) while the 
latter is a non-diffusive momentum flux. The turbulent viscosity is negative with inverse cascade of energy as long as 
rotation is sufficiently strong compared to flow shear (£1 3> A). This is consistent with previous works which showed 
that a turbulent viscosity exists only for highly anisotropic flows [46| or two-dimensional flows [54] . As rotation tends 
to make flow two-dimensional, we expect the turbulent viscosity to be negative. In comparison, turbulent viscosity 
is positive in the opposite limit of weak rotation (£1 A). This is because rotation favors transfer of energy from 
small scales to large scales (inverse cascade) while flow shear efficiently creates small scales via shearing, cascading 
the energy from large to small scales. Even if the eddy viscosity is negative for strong rotation (f2 3> -4), flow shear, 
which transfers energy to small scales, has an interesting effect by slowing down the rate of inverse cascade with the 
value of negative eddy viscosity decreasing as \vt\ oc A~ 2 for strong shear. 

The non-diffusive part of momentum transport (A-effect) can act as a source of large-scale flow, preventing a uniform 
rotation to be solution of the averaged Reynolds equation. A strong anisotropy induced by flow shear Q gives rise to 
non-trivial A-effect even for an isotropic forcing. Note that in the absence of flow shear, the appearance of a A-effect 
requires a source of anisotropy in the system such as an anisotropic forcing in which case the A-effect is proportional 
to the anisotropy in the velocity field [U [33[. Interestingly, our results show that the A-effect scales as A~ 2 \ ln£| 
whereas the anisotropy in the velocity amplitude is given, to leading order, by A~ 4 ^ 3 . Consequently, the A-effect is 
smaller than the anisotropy in the turbulent velocity amplitude. This is because the anisotropy is not simply given 
here but has to be induced self-consistently by the shear during the evolution. In other words, the anisotropy does 
not remain the same at all time, and the resulting A-effect is smaller than the anisotropy in the velocity amplitude in 
the long-time limit. One can also note that the magnitude of the A-effect is not the same in the two cases. 

C. Transport of particles 

The dynamics of particles transport crucially depends on whether rotation is stronger or weaker than flow shear. 
When rotation is stronger than flow shear (f2 3> A), the transport is inhibited by inertial waves, being quenched 
inversely proportional to the rotation rate (i.e. oc i7 _1 ) while in the opposite case where flow shear is stronger than 
rotation, it is reduced by shearing as A^ 1 . It is important to compare this result with turbulence amplitude, which is 
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quenched by shearing even when £1 3> A. This strikingly different behavior between particle transport and turbulence 
amplitude highlights the different roles of waves and flow shear in turbulence regulation; that is, waves mainly affect 
transport by altering phase relation while flow shear quenches both transport and turbulence level, via enhanced 
dissipation. 

Furthermore, in the strong rotation limit (Q 3> A) where the transport of particles is dominated by inertial waves, 
the transport is almost isotropic with only a slight anisotropy - the transport in the direction parallel to the rotation 
is twice larger than the one in the perpendicular direction (see also Eq. (1421) and [34|). However, in the weak rotation 
limit, it is flow shear that quenches particle mixing; the anisotropy in resulting transport can be very large with much 
slower mixing by a factor of A~ 2 / 3 in the direction of shear. The rotation on shear-dominated turbulence weakens 
the anisotropy. 



The calculation of all the turbulent coefficients in the weak shear limit (£ ^> 1) and also of the transport of particles 
in the strong shear limit (£ <C 1) required the evaluation of the integrals of the following type: 



where TUq = (fi ■ k)/fc is the projection of the unit vector in the direction of the wave number on the rotation axis. 
When the domain of integration is unbounded (infinite), the integration over the angular variable of this integral 
becomes proportional to SI -1 , when the rotation rate fl is sufficiently large [see Eq. for details]. This is 

because this integral involves some contribution of order unity (when Jl ■ k = 0) and others of magnitude f2~ 2 . 

However, in realistic situations, the domain of integration in Fourier space is bounded with a minimum wavenumber 
that is permitted in the system (corresponding to a maximum length, for instance the size of the box) in the direction 
of the rotation. If we call this minimum wavenumber k m — min(k x ), we can show that the preceding scaling of fi -1 
is valid only when v 2 k & ^> ^h 2 ^. ^ n the opposite case, the term uJ^ 2 in Eq. ([50]) is always dominant, altering this 
integral to cx fl~ 2 for large rotation rate, with a stronger dependence on Q. 



In this paper, we have performed a thorough investigation of the combined effects of shear and rotation on the 
structure of turbulence, by using a quasi-linear theory. We assumed an external forcing in the Navier-Stokes equation 
which leads to an equilibrium situation where the dissipation (whose effect is enhanced by the shear) is balanced by 
the injection of energy due to forcing. It is useful to recall that there are three (inverse) time-scales in the problem: 
the shearing rate A, the rotation rate £7 and the diffusion rate T> — vk 2 where v is the (molecular) viscosity of the 
fluid and k~ x is a characteristic small-scale of the forcing. The first regime of strong rotation (fi 3> A) has been 
studied in the strong shear (A ^> T>) and weak shear (A <C T>) limits. However, the second regime of weak rotation 
has been considered only in the strong shear (A 3> T>) case, as the effects of both shear and rotation disappear in the 
opposite case. 

While both rotation and (stable) shear flow tend to regulate turbulence, there are important differences in their 
effects, which should be emphasized. Rotation, by exciting inertial waves, tends to reduce turbulence transport more 
heavily than turbulence amplitude while shear flows reduce both of them to a similar degree. That is, rotation (or 
waves) quenches the cross-phase (normalized flux) more than shear flow does [13, Ht1|. Furthermore, in sharp contrast 
to rotation, shear flow induces a strong anisotropic turbulence and transport (e.g. momentum transport, chemical 
mixing, etc.). 

Specifically, in the large rotation limit (|0| 3> A), we have found: 

• The turbulent intensity is reduced only by a strong shear (i.e. in the case of strong rotation and strong shear) 
and in an anisotropic way. 

• As the dissipation decreases (compared to the shear), there is a crossover from a positive to a negative viscosity. 



D. Effect of a bounded domain 




(50) 



V. 



CONCLUSION 



• The transport of particle is reduced by rotation, with a slight anisotropy of a factor 2, largely unaffected by 
shear. 
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In the opposite weak rotation limit (|fi| <C -4), we found that the main reduction is due to the shear with an 
anisotropic turbulence with preferred motion and transport in the plane perpendicular to the shear. Rotation can 
increase or decrease slightly the turbulence intensity and the particle transport, depending on the sign of = fi/A 

Furthermore, we found non-diffusive flux for momentum transport (the so-called A-effect) which transfers energy 
from the fluctuating velocity field to the large-scale flow. In the large rotation limit, this term can appear only for 
an anisotropic forcing. In contrast, in the weak rotation limit, rotation acting together with shear flow was shown to 
give rise to non diffusive fluxes even with an isotropic forcing. 

These results can have significant implications for astrophysical and geophysical systems. For instance, the A-effect 
and/or negative viscosity can provide a mechanism for the generation of ubiquitous large-scale shear flows, which 
are often observed in these objects. Furthermore, the anisotropic mixing of scalars should be taken into account in 
understanding the surface depletion of light elements in stars [551 ] - Finally, we note that numerical confirmation of 
our prediction and the extension of our work to stratified rotating sheared turbulence with/without magnetic fields 
remain challenging important problems, and will be addressed in future publications. 
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APPENDIX A: LINEAR STABILITY ANALYSIS OF THE HOMOGENEOUS SYSTEM 



As Eq. ([2|) is the same as that for a perturbation u about a basic flow Uo, up to the extra forcing term f, our 
study gives some insight into the stability of shear flows in presence of rotation. After summarizing results previously 
obtained by others, we present our results in the cases where the rotation and the shear are perpendicular and parallel, 
respectively. 

The case of the plane shear flow in a rotating frame has been studied by many authors focusing on the stability both 
in the laminar and the turbulent cases. In the case of a rotation vector Jl = £le z perpendicular to the plane of the 
shear flow, [56| proposed an analogy between rotation and stratification [supported by calculation of [57| ] and showed 
that the system was unstable if the vorticity of the shear flow — Ae z is anti-parallel to the rotation and sufficiently 
strong. Precisely, the ratio ft = 2£l/A must lie in the interval [0, 1] for instability. This destabilization of laminar 
shear flow by rotation has a counterpart for turbulent flows where the rotation can stabilize turbulence (by decreasing 
its kinetic energy) or destabilize it, as shown by [58j using a displacement argument. It is interesting to note that 
both Bradshaw and Tritton arguments are pressure- less. However, the Pedley criterion was shown to hold by [59| . 
using stability analysis confirmed by simulations [60] : the cyclonic shear (J7 < 0) is always stabilizing whereas the 
anticyclonic shear (f2 > 0) is destabilizing for weak rotation while stabilizing for high rotation, in agreement with 
Bradshaw criterion. These conclusions are confirmed for a Poiseuille flow, both experimentally [6l| and numerically 
[62| . and for a plane Couette flow (fTl, 63]. The fact that the pressure-less argument gives the exact stability criterion 
is due to the fact that the modes which are dominant in the instability process are naturally unaffected by pressure 
fluctuations . 

The exact solution of the homogeneous part of Eq. (115[) for the velocity v x can be found in terms of generalized 
hyper-geometric function F([ai, ■ ■ ■ ], b 2 , . . . ], x) |65j. Two independent solutions are: 



Xi{t) = F( 
X 2 (t) = tF 
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Here, b — (3 2 £l(£l — 1) is (up to the multiplicative constant f3 2 ) the quantity introduced by [56[ (see discussion in the 
introduction). Figure [3] shows the evolution of these two functions as a function of r. 

Solutions for the other components of the velocity are obtained by using the last two equations of (fTBf : 



Vy = ~- [TX n (T)+p 2 (Cl-l)Y n ( T )] , 

v z = -- [rX n (r) - (Q - l)y„(r)] , 

7 



(A2) 
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FIG. 4: Evolution of the solution Yi (panel [a]) and Y2 (panel [b]) as a function of r for b = —0.5 (circles), 6 = (crosses) and 
b = 0.5 (squares). 



for n — 1 or 2. Here, Y\ and I2 are defined as: 
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The plots of Y\ (r) and F2 (t) are shown in figure 0] 

Figure 0] shows that the eigenfunctions diverge for r — > 00 when 6 < 0. This is because shear flows in presence of 
rotation (perpendicular to the shear flow) is stable only for 6 > 0. This result agrees with [56| and [2{J. We can also 
notice that the solution with 6 > always decays faster than that with 6 < 0. 

In conclusion, we recovered the Bradshaw criterion [561 ]. In our notation, it states that the configuration is unstable 
if B = — 0(1 — Q) < or, equivalently, if Vl = VtJ A lies in the interval [0 , 1]. This result has already been reported by 
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many authors, who showed not only that the maximum destabilization occurs for f2 = 1/2 but also that there is an 
important asymmetry with respect to £l — 1/2 which is not included in the Bradshaw criterion [5^, [66|, [67| • This is 
because Bradshaw criterion can be recovered by a pressure-less analysis: while the pressure does not affect the most 
unstable modes [H| and thus does not alter the instability criterion, is does destroy this symmetry. We can easily 
show that there is indeed asymmetry with respect to ft = 1/2 in our results: even if Eq. (|A1|) for the x-component 
of the velocity is symmetric with respect to O = 1/2 (as it depends only on b = —f3 2 B), Eq. (IA2|) for the other 
components of the velocity are not because of the term proportional to — 1. 



APPENDIX B: DERIVATION OF THE PARTICLE TRANSPORT IN THE LARGE ROTATION LIMIT 

By using Eqs. ([9]), (fT"7|) and (fl9|) . we obtain the transport of particles in the direction of the shear as follows: 

D T = / d 3 k 7 ( 1 + a^F(k) % , (Bl) 

(2lT) i A 2 J UJ 

T (2»)M 2 J 7^1 



where: 

r+oo p -2£[Q(t)-Q(o)] 

C£( k ) = _/ ( 7 + T 2)n/4 exp[^ { a ( O )- g (r)}]dr, (B2) 

/•+oo r p e -2e[Q(T)-Q(o)] 

Mn(k) = y — ( 7 + r 2)«/4 — (t - a) exp[iiv {s(a) - s(t)}] d,T , 

CI = , SI = 3(C) , C£ = , S£ = 3(/<) . 

The expression for is omitted here as it is very similar to that for Dtp* . The asymptotic behavior of integrals 
(|B2|) can be obtained to leading order in Qq 1 as: 

( 7 + a 2)n/4[ 4l/ 2 fc 4 + ^2 ] . ^ 



where a^o = ooqA./ y7 + a 2 - In comparison, the functions /i^ vanish to leading order and are thus omitted here. By 
using Eq. (|B3|) in Eq. (|B1[) . we obtain the following results: 

(B4) 

D T = (2^3 / ( 7 + f 2) 4 ^2 fc 4 + ^2 - 



(2tt) 3 7 v ; ( 7 + a 2 )4^ 2 fc 4 + c^ 2 ■ 

Here, we have discarded all the terms which are odd in a (for example in , the terms proportional to C\ and C3) 
as they vanish after angular integration. After performing this integration, (|B4|) reduces to Eq. |42|) given in the 
main text. 
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